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1 Introduction 

The importance of Sturm-Liouville problems for spectral methods lies in 
the fact that the spectral approximation of the solution of a differential 
equation is usually regarded as a finite expansion of eigenfunctions of a 
suitable Sturm-Liouville problem. The issue of expansion in generalized 
eigenfunctions is a classical one going back at least to Fourier. A relatively 
recent impact is due to the study of wave propagation in random media 
[H [6] , where eigenfunction expansions are an important input in the proof 
of localization. The use of this tool is settled by classical results in the 
Schrdinger operator case. But with the study of operators related with 
classical waves, [II [7], a need for more general results on eigenfunction 
expansion became apparent. Eigenfunction expansions problems for Sturm- 
Liouville problems have been investigated by many authors, see [21 [8] 

In this paper we shall investigate one discontinuous eigenvalue problem 
which consists of Sturm-Liouville equation, 

r(y) := -y"{x, A) + q{x)y{x, A) = \y{x, A) (LI) 



to hold in finite interval (— vr, tt) except at one inner point € (— vr, vr) , 



where discontinuity in u and u' are prescribed by transmission conditions 
ri(y) := aiy'(0-, A) + 022/(0-, A) + a3y'(0+, A) + a4y(0+, A) = 0, (1.2) 

My) := hy'iO-, A) + 622/(0-, A) + bsy'iO+, A) + 642/(0+, A) = 0, (1.3) 
together with the boundary conditions 

^3(2/) '■= cosay{—TT, A) + sinay' {—n, A) = 0, (1-4) 

T^iy) := cos /32/(7r, A) + sin /3y'iTT, A) = 0. (1.5) 

Such problems with point interactions are also studied in [5], etc. Boundary 
value problems which may have discontinuities in the solution or its deriva- 
tive at an interior point are also studied. Conditions are imposed on the left 
and right limits of solutions and their derivatives at an interior point and 
are often called transmission conditions" or interface conditions". These 
problems often arise in varied assortment of physical transfer problems, see 
[His]. Also, some problems with transmission conditions arise in thermal 
conduction problems for a thin laminated plate (i.e., a plate composed by 
materials with different characteristics piled in the thickness). 

It is the purpose of this paper to present a new and somewhat ex- 
tensive family of orthogonal expansions associated with an discontinuous 
Sturm-liouville problem. By using an own technique we introduce a new 
equivalent inner product in the Hilbert space L2(— vr, 0) © L2(0,7r) and a 
linear operator in it such a way that the considered problem can be inter- 
preted as eigenvalue problem for this operator. Using a new approach to 
the now classical theory of Sturm-Liouville eigenfunction expansions, based 
essentially on the method of integral equations. Moreover, we introduce a 
expansion method for solution of carlemans equation. 



2 Preliminary Results 



Let T= 



Denote the determinant of the i-th and j-th 



ai 02 03 04 
61 62 63 64 

columns of the matrix T by pij . Note that throughout this study we shall 
assume that pi2 > and psi > 0. 

In this section we shall define two basic solutions (p{x, A) and x(a;. A) 
by own technique as follows. At first, let us consider solutions of the equa- 
tion (II. ip on the left-hand [— vr, 0) of the considered interval [vr, 0) U (0, vr] 
satisfying the initial conditions 



2/(-7r, A) 



sma, 



dy{-iT,X) 
dx 



cos a 



(2.1) 
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By virtue of well-known existence and uniqueness theorem of ordinary dif- 
ferential equation theory this initial- value problem for each A has a unique 
solution 01 (x, A). Moreover [[8j, Teorem 7] this solution is an entire func- 
tion of A for each fixed x € [— tt, 0) . Using this solutions we can prove 
that the equation (jl.ip on the right-hand interval € (0, vr] of the considered 
interval [— vr, 0) U (0, vr] has the solution u = (p2{x,X) satisfying the initial 
conditions 

y{0, A) = — (p230i(O, A) + p,^^^lp^) (2.2) 

Pl2 OX 

y'{0,X) = ^(pi3</>i(0,A) +pi4^^^). (2.3) 

Pl2 OX 

By applying the method of [4] we can prove that the equation (jl.ip on 
(0,7r] has an unique solution (j)2{x,X) satisfying the conditions (j2.2j) - (l2.3p 
which also is an entire function of the parameter A for each fixed x G (0, tt]. 
Consequently, the function (j){x, A) defined by 

</.(x. A) = {'^i(^''^,\^°^^^ (2.4) 
' ^ <f>2{x,X) forx€(0,7r]. ^ ^ 

satisfies equation (jl.ip . the first boundary condition ()1.4p and the both 
transmission conditions (|1.2p and (11. Sp . Similarly, X2ix,X) be solutions of 
equation (jl.ip on the left-right interval (0, vr] subject to initial conditions 

y(vr,A) = -sin/?, ^^^^ = cos/?. (2.5) 

By virtue of [j8], Teorem 7] each of these solutions are entire functions of A 
for fixed x. By applying the same technique we can prove there is an unique 
solution xi{x,X) of equation (jl.ip the left-hand interval [— 7r,0) satisfying 
the initial condition 

y(0, A) = — (mX2(0, A) + p24^M^), (2.6) 

P34 OX 

y'iO, A) = — (P13X2 + A23^^4^). (2-7) 

P34 ox 

By applying the similar technique as in [4] we can prove that the solu- 
tions Xi{x, A) are also an entire functions of parameter A for each fixed x. 
Consequently, the function x(a;, A) defined as 



xix,X) 



Xiix,X), a;G[-7r,0) 
X2{x,X), a;G(0,7r] 
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satisfies the equation (jl.ip on whole [— 7r,0) U (0,7r], the other boundary 
condition (jl.Sp and the both transmission conditions (|1.2|) and (|1.3|) . In 
the Hilbert Space % = L2[—t^.,^) © L2(0,7r] of two-component vectors we 
define an inner product by 

TT 



<y,z>u--= Pi2 J y{x)z{x)dx + J y{x)z{x)dx 

-TT 

for y = z = z{x) € 7i. We introduce the hnear operator A : 7i ^ V. 

with domain of definition satisfying the following conditions 

i) y and y' are absolutely continuous in each of intervals [— vr, 0) and (0, vr] 
and has a finite limits y(c=F) and y'(c=F) 

ii) Ty{x) G n, riy(x) = T2y{x) = T^yix) = F^x) = 0, 
Obviously D{A) is a linear subset dense in H. We put 

{Ay){x) = Ty{x), x € H 

for y € DiA). Then the problem (jl.ip — (jl.Sp is equivalent to the equation 

Ay = Xy 

in the Hilbert space T-L. Taking in view that the Wronskians W{(j)i,Xux) := 
(pilx, X)Xi{x, X) — (l)'^{x, X)xi{x, X) are independent of variable x we shall 
denote 'Wi{X) = W{(j>i,Xi',x) {i = 1,2). By using ()1.2p and (jl.3p we have 
pi2Wi{X) = P34'W2{X) for each A G C. It is convenient to introduce the 
notation 

w{X) := psiWiiX) = pi2W2{X). (2.8) 
Theorem 2.1. For aUy,z £ D(A), the equality 

< Ay, z >=< y, Az > (2.9) 

holds. 

Proof. Integrating by parts we have for all y,z £ DiA), 


<Ay,z> = pi2 J Ty{x)z{x)dx + pu J Ty{x)z{x)dx 

— TV 

n 



Pi2 / y{x)Tz{x)dx + P34 / y{x)Tz{x)dx + pi2W[y,z;0- 







- Pi2W[y, z; -vr] + P34W[y, z; vr] - p^.^ W[y, z; 0] 

= <y,Az> +pi2W[yo,z;0] - pi2W[y,z; -vr] 

+ P34W[y,z;TT]- puW[y,z;0] (2.10) 
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From the boundary conditions (|1.2p - (|1.3p it is follows obviously that 

W{y,z;-TT) =0 and W{y,z; it) = (2.11) 

The direct calculation gives 

pi2W{y,z;0) = p3iWiy,z;0). (2.12) 

Substituting (12111) and (|2J2]) in ([2T0]l we obtain the equahty (gJl). □ 

Relation (|2.9p shows that the operator A is symmetric and selfa-adjoint. 
Therefore all eigenvalues of the operator A are real and two eigenfunctions 
corresponding to the distinct eigenvalues are orthogonal in the sense of the 
following equality 

n 

Pi2 J y{x)z{x)dx + J y{x)z{x)dx = 0. (2-13) 

-TT 

3 Expansion theorem by Green function's 
method 

To present its explicit form we introduce the Green function 



G{x, s;\) = < 



for — TT < S < X < TT, X, S 7^ 



for — TT < X < S < TT, X, S 7^ 



(3.1) 



for x,s £ [— TT, 0) U (0, vr] where (j){x, A) and x{x, A) are solutions of the 
boundary- value- transmission problem (jl.ip - (jl.5p . It symmetric with re- 
spect to X and s, and real-valued for real A. We show that the function 

TT 

y(x,A) = pi2 j G{x,s;X)f{s)ds + p34 J G{x, s; X)f{s)ds (3.2) 

-IT 

called a resolvent is a solution of the equation 
y" + {X-q{x)}y = f{x), 



(3.3) 
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(where f{x) 7^ is a continuous function), satisfying the boundary-transmission 
conditions (|1.2|) - p.5|) . We can assume that A = is not an eigenvalue. 
Otherwise, we take a fixed number 77, and consider the boundary- value- 
transmission problem 

y"{x, A) + {(A + 7?) - q{x)}y{x, A) = (3.4) 

aiy'(0-, A) + a2y(0-. A) + a3y'(0+, A) + 042/(0+, A) = 0, (3.5) 

6iy'(0-, A) + 62^(0-, A) + 63y'(0+, A) + 642/(0+, A) = 0, (3.6) 

cos a2/(— vr. A) + sin ay' {—it, A) = 0, (3-7) 

cos/3y(7r,A) +sin/32/'(7r,A) = (3.8) 

with the same eigenfunction as for the problem (ll.ip - (11.4p . All the eigen- 
values are shifted through ij to the right. It is evident that r] can be selected 
so that is not an eigenvalue of the new problem. 
Let G{x, s; 0) = G{x, s) then the function 





y{x,X) = pi2 / G(x,s)/(s)ds + />34 / G{x,s)f{s)ds (3.9) 



is a solution of the equation y" — q{x)y = f{x), and satisfies the boundary- 
transmission conditions (jl.2p - (|1.4p . We rewrite (|3.3p in the form 

y" - q{x)y = fix) - Ay (3.10) 

Thus, the homogeneous problem (/(x) = 0) is equivalent to the integral 
equation 

TT 

y{x,X) + X{pi2 J G{x,s)yis)ds + J G{x, s)y{s)ds} = (3.11) 



Denote by Aq, Ai, A2, A„, ... the collection of all the eigenvalues of the 
problem (|l.ip - (jl.4p . and the corresponding normalized eigenfunctions by 

^0,^1,^2, ■■■,^n, ■■■ 

ipn{x)ipn{0 



K{x,0 = Y. 



„ An 

n=0 
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By the asymptotic formulas for the eigenvalues, obtained in previous 
section, the series for H{x,S^) converges absolutely and uniformly; therefore, 
K{x,S^) is continuous. Consider the kernel 



Fix, = + K{x, e) = G(x, + 



n=0 



which is obviously continuous and symmetric. By familiar theorem in the 
theory of integral equations, any symmetric kernel P{x, ^) which is not 
identically zero has at least one eigenfunction [?], i.e., there is a number Aq 
and a function u{x) ^ satisfying the equation 

IT 

u{x, A) + Ao{pi2 J P{x, i)u{i)di + />34 j P{x, i)u{i)di} = 0. (3.12) 

-TT 

Thus, if we show that the kernel has no eigenf unctions, we obtain P{x,^) = 

0, i.e., 

G(x,e) = -V^^^^4^- (3-13) 

n=0 

Hence, to obtain the completeness of eigenfunctions is now easy. It follows 
from the equation ()3.1ip that 

n 

Pl2 I G{X, O^mC + P34 I G{X, O^im = -J^^fnix) (3.14) 

-TT 

therefore, 

TT 

P12 I P{x, CMOd^ + P34 I P{x, CMOdC = (3.15) 

-TT 

1. e., the kernel Q{x,S,) is orthogonal to all eigenfunctions of the boundary- 
value-transmission problem ()l.ip -( fL4]l . Let u{x) be a solution of the inte- 
gral equation (|3.12p . We show that u(x) is orthogonal to all (pn{x). In fact, 
it follows from (f3l2]l 

TT 

= pi2 / u{x)ipn{x)dx + P34 / u(x)(/?„,(x) 
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therefore, 

TT 

= n(x,A) + Ao{/9i2 J P{x,CM0dC + P34 J Pix,CHC)dC} 

-TT 

TT 

= u{x,X) + Xo{pi2 j G(x,e)n(e)(ie + P34 j G{x,i)u{Odi}, 

-TT 

i.e., n(x, A) is an eigenfunction of the boundary- value-transmission problem 
(jl.ip - ()1.4p . Since it is orthogonal to ah (pn{x)-, it is also orthogonal to itself, 
with the consequence that m(x, A) = and = 0. The formula (|3.13p 

is thus proved. 

Theorem 3.1. (Expansion Theorem) If f[x) has a continuous second 
derivative and satisfies the boundary-transmission conditions il.S^) - [L4\ ), 
then f{x) can be expanded into an absolutely and uniformly convergent 
Fourier series of eigenf unctions of the boundary-value-transmission prob- 
lem (E2p-([Il]j on [-1,0) U (0,1], i.e., 

oo 

/(a;) = ^c„99„(x) (3.16) 

n=0 

where Cn are the Fourier coefficients of f given by 

TT 

Cn = pi2 j f {x)^pn{x)dx + j f {x)ipn{x)dx. (3.17) 

-TT 

Proof □ 

Theorem 3.2. (Parseval equality) For any square-integrable function f{x) 
in the interval [—1,0) U (0, 1], the Parseval equality 

P12 / fHx)dx + P34 / fHx)dx = ^ cl (3.18) 
«=o 

holds. 

Proof. □ 
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4 The modified Carleman equation 



We now return to the formula 





y{x,X) = pi2 J Gix,s;\)f{s)ds + J Gix,s;X)f{s)ds (4.1) 





whose right-hand side has been cahed the resolvent. Let 

^ 



y(a;, A) = ^6„(A)(^„(x), c„ = pi2 / /(x)<^„(x)(ix + P34 / /(x)(/?„(x)dx(4.2) 



n=0 

Then, we have 




Cn = P12 j {y(x)" - q{x)y{x}ipn{x)dx + p2,i j {y(x)" - q{x)y{x)}ipn{x)dx 

-TT 

= -A„6„(A) + 6„(A). (4.3) 
Hence, 6„(A) = ^'^'^ ^^'^ expansion of the resolvent is 

y(x,A)=pi2 y G{x,s-\)J{s)ds + p3^ J G(x, s; A)/(s)ds = ^ ^^^4.4) 

-TT 

An important formula can now be derived from the above. Substituting 
the equality 

TT 

Cn = Pl2 J fis)(Pn{s)ds + p34 J f{s)(pn{s)ds (4.5) 

— TT 

on the right-hand side, we see that 

TT 

Pi2 J G{x,s;X)f{s)ds + P34 J G{x,s;X)f{s)ds 

-IT 

OO TT 

= Efr£i'^i2 J fis)Ms)dx + ps, J f{s)Ms)dx}. (4.6) 



n=0 



Since f(s) is arbitrary, 



oo 



G{x,s;t) = 



(4.7) 



n=0 



Thus we obtain 

~ P oo 

Pi2 / G{x,x;t)dt + / G{x,x;t)dx = 




(4.8) 



Put A^(A) = X^o<A„<A number of eigenvalues A„ less than A 
from t (j4.8p the modified Carleman equation 



we get 




dN{X) 
t-X' 



(4.9) 
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